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Abstract 

The decay of the hysteresis loop area of the system, which is obeying a site diluted kinetic Ising model, is 
considered by the disorder parameter using the effective field theory analysis. The exhibition focuses on the 
understanding of external field frequency, amplitude and the site concentration dependency of the hysteresis 
loop area for several powerful treatments. Important characteristics of the hysteretic response, such as fre- 
quency dispersion, effect of domain nucleation phenomenon on the dynamic process etc. has been introduced 
together with well known other characteristics. An attempt has been made to explain the relations between 
the competing time scales (intrinsic microscopic relaxation time of the system and the time period of the 
external oscillatory field) and the shape of the response. As a result of the detailed investigations, existence 
of essentially three, particularly four types of dispersion curves have been propounded. 
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1. Introduction 

When a cooperative many-body system such as a magnet is subject to an alternating external magnetic 
field, the ordinary magnetization, in other words, time dependent magnetization may not respond to the 
external field simultaneously. Hence the system exhibits an ubiquitous phenomenon which is called hystere- 
sis. The phenomenon is a signature of how the system parameters respond dynamically to the external field 
sweep. For the kinetic Ising-type systems, the most familiar hysteresis in the form of Lissajous curve is m{t) 
of the system versus the external magnetic field h(t). Early experimental studies involving ground-breaking 
discoveries in the understanding of hysteresis on real magnetic systems by Steinmetz, dating back to last 
century, there is an empirical law for the hysteresis loops (HLs) areal scaling [![. Also, the hysteresis loop 
area (HLA) has been analyzed as a function of the two open variables of the external magnetic field (ampli- 
tude and frequency) by fitting the data from large-scale (wide-range) to an areal Lorentzian and/or different 
scaling law with from low values of the relevant variables to a monotonically decreasing lawproposed at var- 
ious times in outstanding papers in the past two decades by different researchers [1, S S HI, 0,11, S HI U|- 
Recently, many efforts have been devoted to prediction and experimental verification of the scaling behavior 
of the HLA for thin magnetic films (for a brief review of HLA scaling results see [Hj]). Apart from these, there 
are detailed geometrical description [l| , a statistical theory of the nonlinear relaxation function to describe 
metastable decay 14 1, and hysteresis criterion based on rate competition li|, 16 1. Additionally, the theoret 



ical observations in the relevant references 0, 0, 113] show that the HLs has a nonzero residual loop area in 
the zero- frequency limit (quasistatic limit). The relevant unassailable experimental observation |l8[ is also 
in agreement with the existence of such residuality. Hysteretic behavior is, however, a complicated process 
of a dynamic and nonlinear nature that elude serious treatment, both experimentally and theoretically since 
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the area enclosed by the HL is directly proportional to the energy loss in a magnetization-demagnetization 
cycle. 

HL originates from a competition between time scales of relaxation treatment of the system and the 
periodic external magnetic field. At high temperatures for the high amplitudes of the external magnetic 
field the system follows the external field with a delay which is called phase lag while this is not the case for 
relatively low temperatures and small amplitudes. This spontaneous symmetry loss indicates the presence 
of a dynamical phase transition (DPT) which shows itself in dynamical order parameter (DOP) which is 
defined as the time average of the magnetization over a full period of the oscillatory field (l9| . DPT is 
firstly observed theoretically by using the master equation associated to the Glauber-type stochastic process 
[20| within the framework of mean field approximation (MFA), and it follows that the time dependent 
magnetization satisfies the Suzuki-Kubo mean field equation [2l[ . Tome and Oliveira located a dynamic 
tricritical point (DTCP) where the type of DPT changes across the dynamical phase boundary in their 
pioneering work. However, the first-order dynamical transitions are accepted as artifact by some researchers 
in the subsequent years 22, 23, 24| 25[. Since then, several experimental studies have been put forward to 
explain the dynamic nature of the phase transition and its mechanism. As an example, a DPT occurring 
in high frequency region is represented by Jiang et al. (26| using the surface magneto-optic Keer-effect 
technique for cpitaxially grown ultrathin Co films on a Cu(OOl) surface. The same technique has been used 
by Nowak et al. [27j to explain the temperature and field dependence of the magnetization reversal process 
on Co28Pt72-alloy samples. In addition, as an experimental example of hysteresis mechanism, frequ ency 



dispersion of the HLA for the ferroelectric material Pb(Zro.52Tio.4s)03 has been examined by Liu et al. [11] . 
For a [Co(4A°)/Pt(7A°)] multilayered system with a strong perpendicular anisotropy, DPT is shown by 
Robb et al. [281 ] . In the relevant study, they reported that the experimental nonequilibrium phase diagrams 
are found to strongly resemble the dynamic behavior calculated from simulations of the kinetic Ising model. 
Hence, the strong evidence of consistency between theoretical and experimental studies can be seen in. 
Polyethylene naphthalene nanocomposites have been investigated by Kanuga et al. [23]. After detailed 
analysis, they showed that the material undergoes three critical structural transitions. Up to now from the 
past periods, it can be seen that increasing of the activities focusing on the most current statement on DPT 
by researchers using Monte Carlo (MC) simulations 0, 0, E3, 0, M, HI, 111 , Effective Field Theory (EFT) 
34 , 35 [ and MFA [7, 17, 33|, |36j]. In some of them, several methods were used as comparison. By employing 
the standard Metropolis MC algorithm with periodic boundary conditions, the kinetic pure Ising model in 
a two dimensional square lattice has been simulated by Lo and Pelcovits [30| ], and they found evidence for 
a DPT. After a comprehensive investigation, firstly in the literature showed by Acharyya and Chakrabarti 
that the HLs are asymmetric in the dynamically ordered phases while they are symmetric in the disordered 
phases around the origin . Subsequently, Acharyya 17 1 presented that the dynamic coercive field exhibits 
a power law frequency variation both in the MFA and MC simulation after high precision numerical and 
computational effort. However, the coercive field calculated by using MFA becomes frequency independent in 
the quasistatic limit. This phenomenon is explained by the author constructing the Landau- type double- well 
free energy mechanism with the thermal fluctuations and nucleation process and specified that the coercivity 
can be thought as a limit of metastability. The similar construction has been done between double-peak 
frequency dispersion with the newly defined two parameters which correspond to the domain-nucleation and 
-growth rates respectively by Zhu et al. (33|. In the first one of the two studies published in the recent 
years, Shi et al. [34j | concluded that the system does not contain the coexistence region where dynamically 
ordered and disordered phases coexist. Following the same methodolo gy, a detailed investigation of the 
kinetic Ising system has been done by Deviren et al. in the second one [351 ] . They found that the system 
exhibits the coexistence region which strongl y d epends on the frequency and amplitude of the external 
magnetic field. In the same year, Punya et al. [361 ] represented the phase diagram in (ksT c / J — ho/ J) plane 
for varying frequencies, in order to show the influences of the frequency on the system and they presented 
in depth satisfactory knowledge about the sweeping in time. Additionally, they reported the frequency 
dispersion of the HLA, the remanence and coercivity can be categorized into three distinct types for a 
fixed temperature. However, the article contains no comments about the domain nucleation process. The 
analytically categorization of the HLs propounded by Lyuksyutov et al. [37j |. In addition to this, the two 
other works about the relations between the domain wall motion and the dynamic hysteresis, are published 



2 



by approximately the same collaborators [38L l39j . In the current year, Idigoras et al. [25| have demonstrated 
that the bias field behaves as if it is the conjugate field of the DOP near the DTP for ferromagnets. Also, 
they showed that DOP which has the conjugate bias field, follows a power law behavior at critical frequency 
and its critical exponent is in agreement with the equilibrium MFA value. Despite the fact that the thermal 
and magnetic properties of the quenched disorder systems such as site (or bond) diluted ferromagnets driven 
by a periodically oscillating magnetic field have been examined within the framework of EFT very recently 
in two works 4(J 41 1 but then as far as we know, the maximum lossy behavior of the such systems have 
not yet been investigated. For the site diluted kinetic Ising model, the global phase diagrams including 
the reentrant phase transition are presented by the authors. They showed that the system exhibits a 
DTCP and coexistence region both of which disappear for sufficiently weak dilution [40(. Following the 
same methodology, the relevant collaborators revealed that the impurities in the system has a tendency to 
destruct the first-order transitions and the DTCP for the bond diluted kinetic Ising model [4l|. And also, 
they showed that dynamically ordered phase regions get expanded with decreasing amplitude which is more 
evident at low frequencies. 

It is well known that the EFT is one of the most powerful method that determines the meniscus which 
separates several phases in the relevant planes, based on the use of rigorous correlation identities as a 
starting point and utilizes the differential operator technique firstly developed by Honmura and Kaneyoshi 
|42| . Although the conventional version of the method fails to find an expression for the free energy, due to 
it takes into account the self spin correlations, the method is superior to MFA which neglects the thermal 
fluctuations via neglecting the self spin correlations. Thus, it is expected from EFT to obtain more reasonable 
results than MFA for these systems, as in the case of static Ising model. Therefore in this work, we intend 
to probe the effects of the quenched site dilution process on the dynamic hysteresis characteristics of kinetic 
Ising model in the presence of a time-dependent external oscillatory magnetic field by using the EFT. These 
types of disorder effects constitute an important role in material science, since the quenched disorder effects 
may induce some important macroscopic effect which is still open for inspection. Eventually, for these 
purposes the outline of the article as follows: We briefly describe the formalism and the method used for the 
problem in Sec. 2. Numerical results and discussions are summarized in Sec. 3, and finally Sec. 4 contains 
remarks about our conclusions. 



2. Methodology 

To investigate the dynamical transitions, dynamical symmetry loss (similar to the phrase symmetry 
breaking in the static Ising model), hysteretic response and many other dynamical features, one simple 
choice may be an Ising ferromagnet defined on the two dimensional space which has a coordination number 
z — 3 with a time dependent external oscillatory (in time but uniform over the space) magnetic field studied 
by EFT. For this purpose, we consider the following Hamiltonian, 

H = - J CiC 3 SiSj ~ h(t) 2J CiSi, (1) 

<ij> i 

where J > favors a ferromagnetic exchange interaction of the adjacent sites, Cj is a site occupation variable 
and Si is the spin variable. Site occupation variable can take the values a = which means that the site i 
is empty or has a non-magnetic lattice item, a = 1 if the site i is occupied by a magnetic item and the spin 
variable can take values Sj = ±1. The first summation in Eq. (fTJ is over the nearest neighbor site pairs and 
the second one is over all lattice sites. The Zeeman term describes interaction of the spins with the field of 
the sinusoidal form 

h(t) = h cos(ujt), (2) 

where t is the time and h$ is the amplitude of the oscillatory magnetic field with a frequency w. The system 
is in contact with an isothermal heat bath at given temperature T . So, the dynamical evolution of the 
system may be given by non-equilibrium Glauber dynamics |20j based on a master equation 



T-n((ciSi)) r = -((asi)) T + ((citanh[/3ci(Ej +h(t))])) r , 
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(3) 



where 1/t is the transition per unit time in a Glauber type stochastic process, j3 — 1/ksT and ks represents 
the Boltzmann constant and Ei is the local field acting on the site i and it is given by 



Ei = J^csss. (4) 



5=1 

The inner average brackets in Eq. ([3]) stands for the usual canonical thermal average and the outer one 
(which has an index r) represents the random configurational average which is necessary for including the 
site dilution effects. 

In order to handle the second term on the right hand side of the Eq. ^ one can use the differential 
operator technique 4^, HJ. By using the technique, Eq. Q gets in the form of 

T^{{CiSi)) r = -{{CiSi)) r + ((Cj exp(c i ^iV))) r /(x)| a!= o (5) 

where V = d/dx is one dimensional differential operator and the function f{x) is given by 

f{x) = ta,nh[/3(x + h(t))], (6) 

as for that the effect of the differential operator on a function f(x), 

exp(aV)f(x)\ x=0 = f(x + a)\ x=0 , (7) 

with any real constant a. By taking into account the two possible values of Cj as Cj = 0, 1 we can write the 
exponential term as 

exp(aci) = C{ exp(a) + 1 — Cj. (8) 
By using the expansion Eq. (O in Eq. ((SJ we can obtain 

T^r{{ciSi}} r = -{(ciSi)) r + ((cj exp(-Ej V))) r /(a;)| x= o, (9) 
where the site occupations cf = Cj was used. By inserting Eq. (Q} into Eq. © we get 



T^((c i s. i )) r = -{{ciSi)} r + (^Ci Y\_[c8 exp(JssV) + 1 - c s ]J J f{x)\ x=Q . (10) 



In order to get a polynomial form of the second term on the right hand-side of the equation, by using the van 
der Waerden identity for two-state spin, i.e., exp(6s;) = cosh(&) + Sj sinh(6) where b is any real constant in 
a representative manner, we write the exponential term in terms of the hyperbolic trigonometric functions, 
Eq. (|10p exactly written in terms of multiple spin correlation functions occurring on the right-hand side. 
Thus we get 

T ^r{{ c i s i))r = -{{ciSi)} r + ( ( Cj [cj cosh( JV) + c s s s sinh(JV) + 1 - c s ] \ \ f(x)\ x = - (H) 
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When the product in Eq. (1111) is expanded, the multi site spin correlations appear. In order to make the 
expansion manageable let us handle these correlations with an improved decoupling approximation (DA) 



44| as 

((Cj . . . CjC k S k . . . CiSi)) r = ((Ci)) r . . . ((Cj}} r ((CfeSfc)) r . . . ((cjSjr)) r , (12) 



with 



(d) r = C 

((CiSi)) r = m. (13) 
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In fact, the primitive form of this improvement corresponds essentially to the Zernike approximation |45[ in 
the bulk problem, and has been successfully applied to a great number of magnetic systems including the 
surface problems 0, |46, 47, 48j . Detailed discussion of DA in site dilution problem for the static case can 




be found in [44j . By using the approximation in Eq. (|T2)l and the definitions in Eq. (fl3|) in Eq. (TTT1) we get 

Tin = -m + c(([ccosh(JV) + msinh(JV) + 1 - c\ z )) r j(x)\ x=0 . (14) 

By using the binomial expansion and writing the hyperbolic trigonometric functions in terms of the expo- 
nential functions we get the most compact form of Eq. (IT41 as 

(15) 

'V ~o J 

where 

K =Yj if) Qc p -' y+1 (l - C ) z - p cosh^^(JV)sinh r (JV)/(x)| x= o. (16) 

Eq. ([T5"j) is typical first order ODE but has not an analytical solution because of the right hand-side contains 
transcendental functions. The dynamical equation of motion can be solved by various numerical methods. In 
this work, we prefer to use the fourth order Runge-Kutta method (RK4) to get the evolution of the m(t) by 
regarding Eq. (fT5"|) as an initial value problem. We can mention that the differential equation derived in Eq. 
(|15p extends up to the term m z . Each term in the equation of motion makes contribution to the solution, 
because the value of m, which is calculated at each time step, is iteratively related to the previous m value, 
however the situation is different from the behavior of the equilibrium systems at which high ordered terms 
can be neglected in the neighborhood of phase transition point. 

The system has three dependent Hamiltonian variables, namely frequency of external magnetic field (w), 
amplitude (ho) and the dilution parameter (c). For certain values of these parameters and temperature, RK4 
will give convergency behavior after some iterations i.e. the solutions have property m(t) = m(t + 2tt/lo) for 
arbitrary initial value for the magnetization (m(t = 0)). Each iteration, i.e. the calculation magnetization 
for i + 1 from the the previous magnetization for i, is now performed for these purpose whereby the RK4 
iterative equation is being utilized to determine the magnetization for every i. In order to keep the iteration 
procedure stable in our simulations, we have chosen 10 4 points for each RK4 step. Thus, after obtaining the 
convergent region and some transient steps (which depends on Hamiltonian parameters and the temperature) 
the DOP can be calculated from 

Q = — I m(t)dt (17) 
2vr J 

where m is a stable and periodic function anymore. As cut-off condition for numerical self-consistency, we 
defined a tolerance 

Q\-2^-Q\T ,u <itr 6 (is) 

meaning that the maximum error as difference between the each consecutive iteration should be lower than 
1CP 6 for all step. On the other hand, the HLA corresponding to energy loss due to the hysteresis is defined 
as, 

A = - <j> m(t)dh(t) = h a uj j> m{t) $m(ut)dt. (19) 

There are three possible states for the system, namely ferromagnetic (F), paramagnetic (P) and the 
coexistence phase (F+P). The solution m(t) in convergent region is satisfied by the tolerance 

m(t) = -m(t + w/w) (20) 

in the P phase which is called the symmetric solution. The solution corresponding to P phase follows the 
external magnetic field and oscillates around zero value which means that the DOP is zero. In the F phase, 
the solution does not satisfy Eq. ((20)) and this solution is called as non-symmetric solution which oscillates 
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around a non-zero magnetization value, and does not follow the the external magnetic field i.e. the value of 
Q is different from zero. In these two cases, the observed behavior of magnetization is regardless of the choice 
of initial value of magnetization m(0) whereas the last phase has magnetization solutions symmetric or non- 
symmetric depending on the choice of the initial value of magnetization corresponding to the coexistence 
region where F and P phases overlap. The remanent magnetization (in other words residual magnetization 
which is the magnetization left behind in the system after an external magnetic field is removed) and 
coercive field which means that the intensity of the external magnetic field needed to change the sign of the 
magnetization, can be calculated by benefiting from the HL in order to understand and clarify the behavior 
of the site diluted dynamical system. 

3. Results and Discussion 

The frequency dispersion, the corresponding coercivity and remanent magnetization for various ho/ J at 
T = 0.5T C (0.5T C = 1.045 where T c is the Curie temperature of EFT DA solution with coordination number 
z = 3 in static case) for a pure ferromagnet which has been rigorously depicted in the text, constitutes as 
the starting point for our systematic investigation, are presented in Fig. ([T]). With the guidance of the brief 
report presented by Zhu et al. t 33'|, as a result of evaluations by the conviction that we have reached, the 
HLA frequency dispersion can be classified into (in particular) four types as type I, II. a, II. b and III, and 
plotted with different color in each panel. In the panel (a) of Fig (JXJ) , type I has single peak functions which is 
labeled as peak II (PII) from the asymmetric HLs, II. a where the area suddenly changes from a symmetric to 
an asymmetric HLs and there is pseudo-double peak functions, PI (why it is not taken into consideration for 
this type of curves will be described below) and the PII (which is prominently in the asymmetric region) , II. b 
where there are markedly double peak functions, PI and PII which are corresponding to smooth directional- 
veering frequency point and symmetry loss of the HLs respectively. In the type II. b, as the frequency is 
increased from PI to PII, a dynamic symmetry loss occurs continuously. And finally, we see that the type 
III has a single peak functions PI from the symmetric HLs with a finite value (residual area) in quasistatic 
limit. We also point out that this labelling procedure can be used for all relevant behaviors in different 
parameter dependency, we have also made it for the concentration dependency. Low coercivity is desired in 
the materials, which means low hysteresis loss per cycle of operation. Based on this fact, we have examined 
both coercivity and remanent magnetization. In the panel (b) of Fig. (JT|), type II. a coercivity appears 
only at low frequencies since at a critical frequency (w c ) values, the hysteresis loop is asymmetric and does 
not cross the fixed-axis. In the II. b of the types, again the coercivity appears relatively at low frequencies 
by a difference from type II. a characteristics: The symmetry loss phenomenon takes place in the PII and 
the coercive field has maximal magnitute at that frequency value. In principle, from the relevant veering 
till the symmetry loss the dynamics controlled by so-called reversal field and at the end of this process a 
gibbous-like eminentia appears as seen in the right arm of the mutual coercivities. Type III coercivity is 
small at low frequencies but gets larger and approaches the amplitude at higher frequencies. The type I 
coercivity are not shown, because no coercive field is found due to the asymmetric behavior. In the panel 
(c) for the remanent magnetization which is belonging to the type I collection, we see that there exists only 
an asymmetric loop. For the type II. a and II. b in the panel (c) and (d) respectively, both positive and 
negative remanence appears at low frequencies, but only positive values appear at high frequencies due to 
the symmetry loss. However, they are separated by a difference such as symmetry loss is discontinuous or 
continuous. Type III in the panel (d), both positive and negative remanence values appear at all the shown 
frequencies. 

Shown in Fig. plotted for better visibility, why the type II frequency dispersion curves which are 
categorized by Punya et al. [36(, is classified into two labeled as II. a and II. b as an extra. As can be seen, 
the most remarkable types of hysteresis related to frequency dispersion are these two particular types, which 
contains features from both types I and III. It clearly shows the effect of frequency on the dynamic HLA 
and dynamic phase. Unlike other stochastic methods such as MC simulation technique which takes into 
account thermal fluctuations, there is a very important difference in EFT DA solution of the remanent. 
Given relatively small values of ho, it is possible that there is only PII in the EFT solution. This is because 
of that the equilibrium can be reached by solving the dynamical equation of motion with ho small enough 
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and T < T c , there can be two real stable time series solutions which oscillate around a non-zero value, 
corresponding to two values of stable equilibrium DOP (a positive and negative one). This means that the 
HLs can be asymmetric even in the quasistatic limit. Just a single peak has been encountered along the 
frequency scale, because all the type I curves are members of the purely asymmetric collection. Simply seen 
that the steep jump of the type II. a remanence before the local maxima due to PII. Because of the value of 
mutual remanences almost unchanged up to the symmetry loss, the peak I is not taken into consideration. 
In the type II. b characteristics, while the positive remanence temporally kinks into the deep, negative 
remanence gets into the window slightly at that small frequency region. This phenomenon means that 
continuous symmetry loss exists due to the PII of the HLA dispersion and the eminentia of the coercivity 
at the same time. Similar to the type I, the type III curves which are member of the purely symmetric 
collection, have the same characteristics and behaves in the same way. For a fixed parameter set, the relative 
maximal loop area is observed in this type of frequency dispersion. Also, in the first two figures, the ho 
dependencies of the HL, coercivity and remanence are considered. If we look at both Fig. (JTJ) and Fig. ^ 
as a whole, at first sight we can see that as ui value increases then the symmetry loss phenomenon occurs for 
the specific ho values. This is an expected result, since the increase in frequencies enhances the phase lag 
between time dependent magnetization and the external field signals. Therefore, the asymmetric behavior 
of the HL becomes more prominent, since the time dependent magnetization has less time to follow the 
oscillatory field. Hence, the system can undergo a DPT which requires a small amount of thermal energy. 
This mechanism is due to the competing time scales in such equilibrium systems and also can be seen in 
the global phase diagrams in (ksTc/ J — ho/ J) planes presented by Akinci et al. There is an exact 

consensus in literature that the HL loses its symmetry when the time period of the oscillating external 
perturbation becomes much smaller than the typical relaxation time of the system. Notwithstanding that 
an increasing in ho for a particular oj causes a decrease in sufficient critical temperature T c value for DPT 
along the relevant (/cbT c /J — ho/ J) transition curve, the same operation for a particular T — 0.5T C causes 
an increase in sufficient critical lo c value for the frequency-induced DPT. 

For a better clarification of the relation between the nucleation process and the frequency-induced DPT, 
should be consulted to the hysteretic behavior and its frequency evolution in Fig. ([3]). Hysteresis loops 
in fundamentally different characteristics relating to principal three agents for various ho/ J at T = 0.5T C , 
from the different type of frequency dispersion are shown in the figure. In both panels (a) and (b) shows 
type II. b at ho/ J = 1.2. Type I at ho/ J = 0.3 and III at ho/ J = 3.5 are seen in the panel (c) and 
(d) respectively. Additionally, each panel have an inset which shows magnetization component against the 
instantaneous field component over a time period for selected from the frequency points as experienced 
significant anomalies. As seen in (a) for ho/ J = 1-2 we observe a smooth directional- veering instead of 
symmetry loss at u — 0.164. Most importantly, two branches of the hysteresis curve intersect each other 
at this point and the HLA has the global maxima due to the multiply of its two-axis. Also, the magnetic 
moment remains an approximately constant value. A little bit the left and the right of this point at 
u) = 0.064 and 0.264 respectively, corresponding hysteresis curves are topologically contra-oriented mutually 
and completely different growth mechanism is onset. For the same ho/ J — 1.2 value in (b) as the frequency 
increases starting from the value at the maximum area of the HL the system undergoes a continuous DPT 
due to the incomplete reversal of magnetic moment. For the critical frequency values at w > 0.354 all the 
HLs are in the form of asymmetric. According to times-series of the magnetization and the external magnetic 
field, increasing the field frequency at first, obstructs the saturation of the ordinary magnetization due to 
the decreasing energy coming from the oscillating magnetic field in a half-time period which facilitates the 
late stage domain growth by tending to align the moments in its direction (i.e. the magnetization begins 
to fail following the oscillatory field) and this makes the occurrence of the frequency increasing route to 
DPT at the critical point. Purely asymmetric, left-facing teardroplike HLs which belongs to the type I 
class seen in Fig. ([3]) (c). We observed the maximal area at uj = 0.54 for ho /J = 0.3 and there is no 
DPT due to the increase of the frequency and corresponding relevant time-series of the magnetization which 
oscillate around a non-zero value as seen in the inset. Finally in panel (d) of Fig. ©, can be seen the 
purely symmetric HLs from III and corresponding magnetization oscillation around the non-zero value in 
the relevant inset. As we can see, a lower frequency range, the magnetization can be more saturated and 
this shows itself as the dominant magnetization processes in the different regions of the curve (i.e. reversible 
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boundary displacements, irreversible boundary displacements and magnetization rotation). Furthermore, 
both types I and III show that the HLA increases at a low frequency and then decreases at a higher 
frequency. The most significant distinction of these two types can be seen in the upper and lower frequency 
limits behavior, considering the differences between types I and III at oj — >• 0, for type I the HLA goes to 
zero, while for type III goes to a finite non-zero value. There are no agents from the type II. a, just because 
its representative HLs's frequency dependent evolution is similar to the panel (b) by a difference that there 
is no directional- veering scenario while as the symmetry loss occurs discontinuously. 

As can be understood from the results described above, two peaks observed PI and PII which corresponds 
to the resonance of symmetric and asymmetric HLs with the period of external field, respectively by means 
of the EFT calculations. According to Gunton and co-workers [49(, the dynamics may consist of domain 
nucleation and/or domain growth such systems that we examined in this study. From this point of view, 
in the symmetrical region of the HLs, both the initial domain nucleation and the late stage domain growth 
are onset, hence the system dynamics is controlled by the combined domain nucleation-and-growth process. 
Esoterically, the type II. a and II. b mechanisms are different from each other. Since for the asymmetric 
HLs region the DOP can be always well above or below zero, at any time during the system evolution we 
can observe most magnetic moments having the same direction. The late stage domain growth is relatively 
inhibited, so the dynamics may be mainly domain nucleation. Lyuksyutov and co-workers explained in detail 
these process mcntioned-above by introducing four characteristic fields [37J . According to their declaration, 
the process of magnetization reversal can be divided into two stages: Nucleation of domains with opposite 
magnetization and growth of these domains. Which process dominates depends on system parameters 
and their history. Our observations are in accordance with these explanations and arguments with some 
unpretentious improvements. 

As one of the most striking aspects of our work is that the dispersion curves evolve from the type I to 
III hierarchically as site concentration decreases for every ho/ J values. As a beginning, we need to dart 
a glance at the relevant phase diagrams in Fig. (j4|) which is plotted for the principal three agents for 
various values of the amplitude ho/ J of the external oscillatory field and the site concentration c so as to 
examine the evolution and underlying mechanism. The agents which could respond to the site dilution the 
most quickly are selected. These are (a) type I at ho/ J — 0.3, (b) II. a at ho/ J — 0.52 and (c) II. b at 
ho/ J — 1.2. The diagram for the type III agent is not shown, because no critical temperature is found due 
to the asymmetric behavior. A portion of the curves along the selected route at T = 1.045 are remains 
the purely asymmetric or remains the purely symmetric region and a small amount of them are signify 
symmetry loss. In order to focus on this interesting disorder-induced evolution we depicted the frequency 
dispersions for ho/ J = 0.3,0.52,1.2,3.5 with the concentration values c = 1.0,0.9,0.8,0.7,0.6,0.5 in Fig. 
([5]). As seen in the figure, at T = 1.045 for a particular disorder parameter, although, a weak increasing 
for ho/ J causes an increase in u) c , & decreasing of the site concentration causes an increase in lo c and vice 
versa. This is an expected results, since as the lattice sites are diluted more and more then the dipole-dipole 
interaction-induced energy contribution gets smaller. The quantity of the energy contribution which comes 
from the dipole-dipole interaction corresponds to the self- volition to remain in the dynamically symmetric 
phase of the system. Hence, the system can stay in the paramagnetic phase in the smaller frequency 
until the frequency of the external field becomes greater than the relaxation time of the system. In order 
to probe the concentration dependence of the dynamic nature of the system, we plotted the HLA versus 
dilution parameter c in Fig. ([6]) for various frequency values (a) ui = 0.01, (b) 0.1, (c) 0.5, and (d) 1.0 
at T = 1.045. Correspondingly, the unmediated concentration dependency of the coercive field and the 
remanence magnetization plotted in Fig. ([7]) and © respectively for the same parameter set acceptances. 
In accordance with our anticipation, as the concentration c of magnetic atoms decreases then the system puts 
under the command of the symmetric behavior. In other words, decreasing of the concentration facilitates 
the symmetry tendency of the system. As a new characteristic behavior, such as discontinuous symmetry 
losses are also available in this anthologia for ho — 0.3 at uj = 0.01. We can also observe the existence of the 
continuous symmetry loss and correspondingly above-cited peak I and II (related eminentia). As a general 
trend, coercivity is small at low site concentrations, but gets larger and approaches the amplitude at higher 
concentration. As seen in Fig. ([7]), increasing the frequency causes the behavior gradually becomes linear and 
in Fig. [8]both positive and negative remanence appears at low frequencies, but only positive values appear at 
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high frequencies due to the symmetry loss. The discontinuous symmetry loss, also the kinks corresponding 
to the relevant eminentia etc. can be clearly seen. Moreover, the low amplitudes concentration dependencies 
in Fig. ©, Q) and flH}, which appear at small w values such asw = 0.01, exhibit a discontinuous symmetry 
loss and the new type II. b curves reproduce hierarchically from II. a curves as increasing frequency. This 
phenomena can seen as the evolution from the sudden symmetry loss to coercivity which has two peaks in 
Fig. (JTJ) and the steep jump behavior to the kink of the remanence in Fig. ©. All statements concerning 
the mechanism of this behavior, also apply here. Furthermore, HLA goes to a finite value in the limit of 
c — > 0. It means that, no matter the system exhibit a long range ordered phase or not, for all concentration 
values we can observe the HLA, correspondingly coercivity and remanence for concentration values even 
below its critical value c* (namely site percolation threshold) . 

Concretely speaking, in order to determine the type of the frequency or disorder-induced symmetry loss 
(continuous or discontinuous) , we appeal to hysteresis for its topological evolution in respect to the relevant 
parameter for all of the investigations mentioned-above. If the evolution exhibits the directional-veering 
phenomenon this type of the loss is classified as continuous otherwise it is assumed to be discontinuous. 
There exists also, reasonable grounds of this classification such as extremely different nucleation mechanisms 
underlying this two different transition types. 

4. Conclusion 

In this study, we have investigated the hysteretic behavior and dynamic nature of the critical phenomena 
which is observed for a site diluted Ising ferromagnet which has been defined with coordination number 
z = 3 driven by an external oscillatory magnetic field by means of effective field theory based on a standard 
decoupling approximation and the time evolution of the system has been presented by utilizing a Glauber 
type stochastic process. Our starting point for the systematic review was to examine the frequency disper- 
sion, coercivity and remanence coordinately for various field amplitude values at finite temperature below 
the Curie point of the static case on a pure Ising ferromagnet. In addition to the classification set forth in 
the literature, we categorized the frequency dispersion into four different types and labeled them according 
to the classification procedure which has been given in the text. We supported the relevant classification 
with hysteresis characteristics as well as the magnetization time series. 

After that, in order to examine the effect of the dilution of site concentration on the dynamic nature 
of the system, firstly we presented the phase diagrams for the selected agents from three distinct types 
for various concentration values c in the (fcgT c / J — lu) planes. In the immediate aftermath, the disorder 
induced evolution of the agents have been investigated at finite temperature corresponds along the route at 
T = 1.045 in the phase diagrams. According to our findings, a weak decreasing of the site concentration 
causes an increase in uj c and vice versa. For a better view, we have depicted disorder dependency of the 
HLA, coercivity and the remanence for the aforementioned agents for various frequency values at T = 1.045. 
Also, our observations were the same, i.e., the relevant different type of dependencies has been observed 
similar to the frequency dependency classes. 

EFT takes the standard MFA predictions one step forward by taking into account the single spin corre- 
lations which means that the thermal fluctuations are partially considered. Although all of the observations 
reported in this work shows that EFT can be successfully applied to such nonequilibrium systems in the 
presence of quenched site disorder, the true nature of the physical facts underlying the observations dis- 
played in the system (especially the origin of the different distinct types of the parametric dependencies of 
the HLA related to the nucleation process and the HL residual area in the quasistatic limit) may be further 
understood with an improved version of the present EFT formalism which can be achieved by attempting 
to consider the multi spin correlations which originate when expanding the spin identities. We believe that 
this attempt could provide a treatment beyond the present approximation. 

In conclusion, we hope that the results obtained in this work would shed light on the further investigations 
of the dynamic nature of the critical phenomena in disordered systems (e.g. scaling of the HLA by using 
the wide-range Lorentzian scaling function with four variabless in the effective-field case) and would be 
beneficial from both theoretical and experimental points of view. 
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Figure 1: The frequency dispersion of (a) HLA, (b) coercive field and (c)(d) the remanent magnetization for various ho/ J 
at T = 0.5T C , which are classified into four types of formal care: Type I at ho/ J = 0.15,0.2,0.25,0.3, II. a at ho/ J = 
0.52,0.54,0.56,0.58, Il.b at h /J = 0.7,0.8,0.95,1.2 and type III at ho/ J = 2.0,2.5,3.0,3.5 respectively. The lines serve as a 
viewing guide. 




Figure 2: (a) Type I, II. a and (b) Il.b, III remanences at T = 0.5T C zoom in as plotted again in terms of visual clarity, 
numbers accompanying each curve denote the value of the amplitude of the external oscillatory field. 
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Figure 3: The HLs in different characteristics for various ho/ J at T = 0.5T C which are located along the frequency dispersion 
representatives from (a)(b) type II. b, (c) I and (d) III are shown in the figure. There is no HLs for the representative which is 
selected from type II. a, because its behavior is exhaustively depicted in the text. Each insets show magnetization component 
against the instantaneous field component over a time period for the related frequency values. 




Figure 4: Phase diagrams for the lattice which has been defined with coordination number z = 3 in (kgT c /J — oj) plane for 
the principal three agents stated in the text. The numbers accompanying each represent value of the site concentration. 
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Figure 5: Disorder-induced evolution of the (a) type I, (b) type II. a, (c) type II. b and (d) type III representatives frequency 
dispersion at T = 1.045. The numbers accompanying each represent value of the site concentration. 
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Figure 7: Variation of the coercivity with site concentration c for the selected agents from four different types of the frequency 
dispersion curves at constant temperature T = 1.045, for the values of the frequency ui = 0.01,0.1,0.5 and 1.0. The numbers 
accompanying each curve denote the value of the amplitude of the external oscillatory field. 
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Figure 8: Variation of the remanence with site concentration c for the selected agents from four different types of the frequency 
dispersion curves at constant temperature T = 1.045, for the values of the frequency uj = 0.01,0.1,0.5 and 1.0. The numbers 
accompanying each curve denote the value of the amplitude of the external oscillatory field. 
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